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Abstract 

We solve a radial Schrodinger equation for the case of a multichannel square well 
plus an exponential potential in one of the channels. The solution is obtained by 
summing exactly the infinite terms of the perturbative series for the evolution operator 
of the system. Analytical expressions for the scattering matrix and the scattering 
length are given, a Deser-like formula which connects approximate ground energies 
and scattering length is obtained for a particular case. Possible physical applications 
of these results are discussed briefly. 
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1 Introduction 



The interactions of pions with nuclei or with other pions have received considerable at- 
tention since the sixties as possible tools for the study of various aspects of the nuclear 
structure or in order to check the validity of diverse effective theories. 

In any process involving bound pions, it is crucial to know the value of the pion atomic 
wave function in the nucleus. The pions are generally trapped in atomic orbits with large 
radial quantum numbers, from which they are very unlikely to be absorbed by the nucleus. 
However, it is possible they cascade downwards, first making Auger transitions between 
higher levels and subsequently by X-ray transitions. It has been known since long that 
for nuclei with Z < 10 some of the pions succeed in reaching the Is orbit before being 
absorbed by the nucleus. The absorption rate is roughly proportional to the probability of 
finding the meson inside the nucleus. This probability for the Is state pions can be altered 
significantly both by the pion-nucleus strong interaction and by the finite charge radius 
of the nucleus. Early numerical studies and experimental data in pionic 2p-ls transitions 
(for example ||l], ||, ^ ) made clear that standard perturbation theory was not adequate for 
the calculation of nuclear pionic wave functions, energy level shifts and widths, even for 
relatively light nuclei. 

Two experimental effects relevant to the study of the interactions of pions with nuclei 
have recently appeared in nuclear physics: the existence of deeply bound narrow pionic 
states in large Z nuclei and the nuclear halo. 

The deeply bound pionic states in Pb, recently discovered at GSI (Q) and previously 
predicted in ( ^ ) , have stimulated a renewed interest in the nature of the s-wave part of 
the pion-nucleus optical potential. For such states to exist with sufficiently long lifetime, 
there must be a subtle cancellation between the vr" coulomb and strong interactions in 
the bulk nucleus. The effective s-wave vr-nucleon repulsion must roughly compensate 
the attractive Coulomb force in the interior of the nucleus. In this way the Is pionic 
bound states in heavy systems may provide strong quantitative constrains on the s-wave 
TT interactions with nucleons. 

The spectra of deeply-bound pionic atoms have been calculated by finding the numer- 
ical solution of a Klein-Gordon equation with a standard pion-nucleus optical potential 
(Q). The widths of the deeply bound states were found to be narrow even for the Is 
state of very large Z. As an alternative to Auger transitions, "pion transfer" reactions 
were proposed by the same authors as a practical way of producing deeply bound pionic 
states. An incoming sufficiently energetic neutron incident on a target nucleus changes 
into an outgoing proton emitting a vr", which then gets captured by the target nucleus 
and forms a deeply-bound pionic state. The process n pTr~ does not take place in free 
space. Hence, at least one of the three particles must be necessarily off-shell, i.e. the vr" 
with highest probability. The {n,p) differential cross section can be written, under the 
Born approximation and supposing plane wave functions for the incoming and outgoing 
neutron and proton, as 
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where ^{q) is the Fourier transform of the pionic wave function. It is then important to 
have rehable methods for the computation of this wave function. 

The neutron halo is a relatively recent discovery in nuclear physics, ([0, §]) which affects 
neutron rich nuclei. This halo effect, which can maintain a significant neutron density even 
at distances of several nuclear radius, can be relevant in diverse neutron-capture reactions 
by neutron-rich nuclei. The neutron halo, being essentially repulsive for vr" and existing 
over a large distance, could favor the existence of exceptionally narrow pion bound states 
for relatively low Z nuclei or at least modify substantially the properties of such pionic 
atoms. Although numerical calculations of Schrodinger or Klein Gordon equations with 
arbitrary complex potentials can be performed easily, it is also convenient to have exact 
analytical solutions, even for simplified models, in order to study systematically such 
effects. These solutions can in fact help us to understand some general properties and 
serve as a benchmark for perturbative approximations and numerical calculations. The 
characteristics of the pion transfer reactions also make highly desirable to count on a 
multichannel framework as an alternative to optical potential methods. 

In this work we present the technical details of the solution of a Schrodinger equation 
with a multichannel well potential plus a Bethe potential in one of the channels. In future 
publications (||9|) we intend to use such a potential, complemented with a Coulombic 
potential in an external region, to describe the general properties of pionium and pionic 
bound states in exotic atoms as a function of A and Z: pionic wave function, level shifts 
and widths. The exponential Bethe potential can be used at will to parametrize locally the 
quadratic electromagnetic potential inside an uniformly charged nucleus of arbitrary size 
or the decaying potential originated by an extended neutron density. The shortcomings of 
using a non-relativistic formalism are compensated by the advantages of having an exact 
analytical multichannel computation. 

Bethe was the first (|]10|) to use an exponential potential to make quantitative calcu- 
lations and explain various general features of cross sections for processes involving nuclei 
and nucleons. He used a function of the form V{r) ~ to describe the complex absorp- 
tive part of the short range nuclear potential. This choice was motivated essentially by the 
possibility of solving exactly the corresponding Schrodinger equation. The same potential 
has also been employed (|11]) to estimate nuclear transmission functions. The model com- 
prising two square-well potentials coupled by a square-well interaction has been used by 
many authors (see for example |12, 13|) in simple specific models of nuclear, atomic and 
molecular physics or as a test of various reaction theories ( ) . By comparing the results 
obtained from reaction theory approximate calculations with the exact analytic solutions 
one gets information about the rates of convergence of the various formulations. 

The work presented here combines multichannel square-wells with the exponential 
Bethe potential. The method which has been used, i.e., summing over the infinite terms 
of a perturbative series for the evolution operator of the system, has at least the following 
advantages. First of all it is of intrinsic interest; it is, to our knowledge, the first time 
that such a method has been used for the solution of a second order Schrodinger equation. 
Moreover it is easily generalizable to other multichannel systems and allows the straight- 
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forward introduction of constant non-local derivative potentials. Finally the evolution 
operator formalism is specially suited, by its own nature, for systems where internal and 
external potentials coexist, as well as the use in conjunction with R-matrix methods. 



2 Solving the Schrodinger equation. One channel case. 

Consider the radial s-wave Schrodinger equation for a central potential of exponential 
profile: 

^ + (q^ + pe-'^)y{r)=0. (1) 

The constants q, p and k are defined by the expressions q'^ = k'^ — 2mV]\f = 2m{E — Vn) 
and p = 2mpQ. The two quantities m and E are respectively the mass and the energy of 
the system; is a constant term which is introduced for convenience and y{r) = rR(r), 
where R{r) is the radial part of the wave function of the system. 

One can reduce Eq.(||) to a first order system of two equations, by defining a two 
dimensional vector ^(r) ( formed by the function y{r) and its derivative) as follows: 



y{r) 

With this definition, Eq. ([l|) is equivalent to the system: 



dr 

where: 



(Ho - pe-^^y) ^{r) , (3) 



We perform a linear transformation ^ = u^d in order to diagonalize Hq. We arrive 
at the equation: 

i-T^ = H{r)^D{r) = {Hd - PDe-^'^W)^D , (5) 
ar 

where pD = pf^Q and 

W = 2iq u-^Vu = }^ ; Hd = i u-^H^u = (^~^ ■ (6) 



iq —iq 
1 1 



(7) 



Note that in the process of passing from the second order equation to the first order 
one, the Hamiltonian has become " squared rooted" : the eigenvalues of Hq are the roots 
of the "operator" q^ and the full Hamiltonian fullfills the relation: 

Hd - PDe-^''W\^ = q^ + pe-^'- (8) 
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We have reduced our original problem to the resolution of a "time-dependent" Schrodinger 
equation with a non hermitic Hamiltonian. The evolution operator U of the differential 
system is given by the path-ordered exponential: 



C/(r, ro) = Pexp 



dsH{s) 



ro 



(9) 



At least formally, it is always possible to give a perturbative expansion for the operator 



U: 



i7(r,ro) = C/(°)(r,ro) + ^C/(")(r,ro); U'-^\r,ro) = e^p[-iHD{r - ro)] . (10) 

n=l 

The n-term U^"^^ in the expansion is given by the integral 

U^^\r,ro) = {-iT j^dTn. . .dnU\r,Tn)W{Tn) . . .U\T2,n)W{Ti)U\n,ro) . (11) 

In the previous formula we used the notation: W{t) = — p^i exp(— At)VF . The 
domain of integration is defined by: 

r = r > r„ > . . . > ri > ro. 

One can show (following the same arguments as in [jl5|) that the computation of U 
in the limit when r — > oo suffices in order to obtain the general evolution operator. The 
evolution operator for a finite interval (r, ro) can be expressed as the product of two 
operators defined respectively at the points r and ro, which are essentially instances of the 
asymptotic U. 

We get the following expression for the elements of C/*^"-' in a basis of eigenvectors of 
Hd through elementary manipulations of Eq.([ll|): 

= {-iT exp {-i{E,r - EaVo)) x 

X ^ / (i"Texp(z(rna;bfci + • • • + Ti6Jfc(„_i)a))VF;,fci(r„) . . . 14^fc(„_i)„(ri) (12) 

fcl,...,fc(n-l) ^ 

where Wfcifc2 = -E-fci ~ -E'fc2 i E^ denotes one of the two eigenvalues of Hq. The matrix W 
is such that W"^ = 0. It has a dyadic form iWij = AiBj) and hence the product of two 
contiguous factors in the Integral (|T^) is 

W^jW.i = WuWjj. 

By applying repeatedly the previous relation, Eq.(|l2|) becomes {Bf^ = Wkk) '■ 

Ula^ = i-^T exp(-i(i?fer - Earo))Wba xY^Bkl--- i?fe(„-l)^fcl...fc(n-l); (13) 

-4fci...fc(„_i) = / d"Texp(r„(iwfefci - A) . . . -|-ri(iu;fc(„_i)„ - A))(-/3d)"- (14) 
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The multiple sum runs independently over all the eigenvalues of the Hamiltonian. The 
limit r — > oo is understood in the domain of integration. One can use then the following 
identity 

In{uJi, . . . ,ujn) = / ■■■/ d^xexpC^uJiXi) 

(-irexp(roEiWi) 



<^n{u^n + t^n-l) ■ ■ ■ {uJn + + • • • Wl) 

in order to arrive at the expression: 

^fci...fc(n-i) = — ■■ 7 — T- T- 7 TTTT- (-'^^) 

Inserting the expression (^) in Eq.(|l^), we obtain for the n-order matrix elements of 



U: 



(n) ^ ^_iE,{r-ro)y^^^^^_i ^^n^-nroX ^ _ ^ 

IWba - nX 



{iujbki - A) . . . {iujbk{n-i) - {n- 1)A) 



e 



iiOba - nX ^^bk - mX 

if.}, — T) \ J- J- 



iujba - nX .^j-^ (-iwfti + mX){-iwb2 + mX) 

. ,-..,.-..1^ ^'"°^)^'"°)" 'Mxr' , (J,) 

-iwfca/A + n (-iwfci/A + l)(„_i)(-zu;62/A + l)(„-i) 

In the third line, we have used the fact that Bi + B2 = and that the expression 
BiWb2+B2Wbi = W12 holds for any value of b. The expressions (a)(b) = a(a+l) . . . (a+n— 1) 
are the Pochammer symbols. The matrix elements of U are given by a convergent series 
which can be summed readily into generalized hyper geometric functions. We have for 
example for f/n, C/12: 

^\ X^ ^ ) n\{-^W^2/X\^) 



00 / 

U^2 = e-^^('--o)^rP£!^ 



-At'o 



n=0 



V A2 / -rai2/A + n (n - 1)! (-iit;i2/A)(„) 



e-^^('--o) (Un - oF, (l - ^e-'A) . (18) 
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The expressions for the other two elements of the matrix U are similar. The matrix U in 
the limit r — > oo can be written as: 

U{r^oo,ro) = exp{-iHDr) Urediro), (19) 

Ured{ro) = exp{iHDro) Us{ro) , (20) 

where the matrix Us is: 

TT _ f oFi{a,z) oFi{a,z) -oFi{l + a,z)\ , . 

{oFi{-a,z)-oFi{l-a,z) oFi{-a,z) J' ^^'> 



The variables appearing in Eq.(^l|) are: a = —iwi2/X, z = PDW12/ X"^ e^p{—Xro). More 
explicitly: a = 2iq/X, z = —p/X'^ exp(— Aro). 

The trace of the Hamiltonian in Eq.(^) is zero. This property alone guarantees that, 
although U is non unitary, its determinant det Us = 1- Hence a relation for the absolute 
values of the hypergeometric functions appearing in Eq.(|2l|) can be obtained immediately 
by-product. 

For finite r the evolution operator for Eq.(^) is given by ( [1l5| ) : 



U{r,ro) = Ured{r)-'Ured{ro)- (22) 

The evolution operator for Eq.(|3|) is obtained by applying an inverse transformation 
with the matrix u defined by Eq.(|^: 

U{r, ro) = u Ured{r)-^Ured{ro) u~\ (23) 

The general solution of Eq.(|l|) is given, in terms of initial conditions at ro, by 

y{r) = U2i{r,ro)y'{ro) + U22{r,ro)y{ro). (24) 

Thus, a regular particular solution, that is, one which vanishes at origin, is given by 

y{r) = C U2i{r,0) , 

where C is an arbitrary constant. 

We have supposed throughout the computations that the parameter A is bigger than 
zero. However the solution is still valid for any value of A by analytic continuation . One 
can check that the function (p4|) is effectively a solution of our equation and coincides with 
the solution given by Bethe (|l(]f|). 

2.0.1 The scattering matrix and scattering length. 

The solution to the radial Schrodinger equation may be written as 

yir)=y+{r)+Siq)y'{r) , (25) 
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where y~{r) and y^{r) denote appropriate ingoing and outgoing solutions to the Schrodinger 
equation with suitable asymptotic behavior at infinity. All the relevant information about 
the system as bound states and scattering parameters can be obtained from the function 
S{q). One can give an explicit formula for S{q) for the case in which the exponential 
potential extends from zero to infinity. In this case the asymptotic states can be read 
directly from Eq.(|l9|) where the matrix Ho is diagonal (Eq.(|6|)). We can write, taking 
into account that ^ = u^o, 

y{r) = + ^'D,2(r) . 

In order to have a function y(r) regular at the origin, the functions ^D,i{f) and ^d.2{^) 
have to be chosen as 

It follows from Eq.(p!9|) that the asymptotic behavior of these functions is 

^D,2(r ^ cx)) ~ Ae-''^''; ^d,i ~ Be''^''. (27) 

The functions ^di,2 are called the "Jost" functions of the system. The asymptotic behav- 
ior of y{r) is: 

y{r) = fe*^'' {Ured,im - Ured,12m + e"'^" (UredM^) " UredM^))] ■ (28) 

Ziq L J 



Given Eq.(|25|) and the expressions for the elements of Ured in Eqs.(2C,21), we obtain 
the following expression for the scattering function: 

5(,) = -£^ll-|?§^. (29) 

oFi(l + 2iq/X, -p/y) 

The bound states correspond to poles of S{q) and the eigenvalue equation becomes 

oFi(l + 2ig/A,-p/A2) = . (30) 
The scattering length A is defined by the limit 

— - = — lim - tan (arg Siq)l2) . 
A q^oq ^ ' 



We can compute explicitly A in this case and from Eqs.(p9|) we deduce that: 

1^.,. 1 oFi(l + 2ig/A, -p/A^) - oFi(l - 2ig/A, -p/A^) 

A \''%q^F^{l + 2iq/\-p/\'^)+oFi{l-2iq/\-p/X^)- ^ > 

The hypergeometric function QFi{a, z) is analytic both in z and in a. We can evaluate 
readily the limit in terms of its derivative with respect to the parameter a: 

1 2 5aoFl(l,-p/A2) 



A A oi^i(l,-p/A2) 



(32) 
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The value of the derivative daoFi(a, z) for a = 1 can be computed deriving term by 
term the series of the hypergeometric function and inserting that value; the result is: 



9„oFi(l,^) = -E^E— (33) 
, nW. '-^^ ni 

n=l m=l 

This series can be summed in terms of Bessel functions: 

daoFi{l,z) = Jo {2i^z) {logi^z + C)- ^No{2i^z) ; (34) 

Jo{z) and Nq(z) are the Bessel functions of the first and second kind respectively. In 
Eq.(^) C is the Catalan number (see for definitions). 



It follows from the eigenvalue equation (pO) and from the formula (32) that we can 



obtain a Deser-like formula |17] which connects the approximate first eigenvalue shift in 
TT-mesic atoms to vr — scattering data. If qq is a sufficiently small root of Eq.(|30|) then 
one finds approximately: 

oFi(l,z) + ^5aoFi(l,z) ~0 (35) 

or 

2 daoFi{l,z) 

Compare this expression with the Deser formula which, under rather general assump- 
tions for the form of the nuclear interaction, can be written as 

— oc ^0 (37) 

q 

where, q and Aq are respectively the energy of the s-state Coulomb bound states and the 
shift produced by the short range nuclear interaction, Aq is the scattering length defined 
in the absence of the long-range (Coulomb) field. This formula can be obtained treating 
the nuclear interaction in first Born approximation. Sophisticated extensions of the Deser 
formula and of Eq.(^) can be obtained by exploiting the analytical properties of the 
scattering function and the solutions of the Schroedinger equation; further development 
will appear in a independent publication. 



3 The multichannel case. 

We proceed now to the generalization of Eq.(|l|) to a multichannel case. Although for 
simplicity and economy of notation we will only work out the two channel case, the ge- 
neral case is obviously analogous. We consider the following matrix equation for the two 
component vector y(r): 
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(2) (2) 

where Hq and are the matrices: 



r{2) _ n_/TA. _ T/(2)=n 



H',> = 2m{VN-Ey, V^o'^^^I^Q q j • (39) 

In Eq.(|39|) m and E are this time the diagonal matrices: m = Diag{mi,m2), E = 
Diag{Ei, E2); Vat is an arbitrary constant matrix. 

Eq.(|3^) can be used to describe a two-channel system where only one of the channels 
is "charged" in the sense that it feels the exponential potential. In realistic applications 
the free parameters of such a potential, p and A, can be adjusted to approximate, at least 
qualitatively, the desired potential. For more quantitative descriptions, the full range of the 
potential can be divided into a number of subintervals and the potential itself approximated 
in each of them by the exponential function. In this way an arbitrary degree of accuracy 
can be obtained. The evolution operator formalism, due to its factorizing property, is 
especially well suited for this numerical approach based on the division into intervals. 

As we did with the one channel case, the second order system given by Eq.(|38[) will 
be transformed into a four dimensional, first order system, whose evolution operator will 
be computed. It is convenient to perform first a diagonalization of Hq ' by applying a 
transformation u: 

H'g> =u-^H^^\ = Diag{Hi,H2); V^^^ = u-^V^''\- y{r)=uyi{r). (40) 

Defining a vector ^'(r) containing y/(r) and its derivative, one arrives at the four- 
dimensional, first order system: 

^-(^«---^'^)*^ *M.(*X)^ (41) 
J^g'V ,,,^0 



One proceeds then to a new diagonalization of the 4x4 matrix Hq. 

Hd = w~^How] W = w'^Vw\ ^' = w^d. (43) 
We are finally left with the system 

dr^D = (Hd - pe-^'w) ■ (44) 

The matrix Ho is 

Hd = Diag[ki,k2, -ki, -k2], ki^2 = \[Hi^2- (45) 
The explicit form of the matrix w is 
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where q is the following matrix 



The 4x4 matrix W is the tensor product of two orthogonal vectors W = a f^i b and 
a-b = 0. The elements of each of these vectors are essentially combinations of the elements 
of the matrix u: 

w _ 1 f —U21 U22 U21 U22 

\/2 det u V ' fei ' A;2 ' h 



V2 det ^ 



(ui2,uii,ui2,nii) . (48) 



Some essential properties of W, that will be used later, are: = 0, Tr{W) = 0. The 
diagonal elements of W are given explicitly by 

(Wbb) = 1/2 {C3/k2,Ci/h,-C3/k2,-Ci/h), 

with 

Ci = U11U22/ det n; C3 = -U21U12/ det u; Ci + C3 = 1. 



The computation of the evolution operator for the two channel version of Eq.(44) 
follows the same steps as the one channel case which we have seen in the previous section. 
In this case, because of the conventions we have used, the factor (— i)"" appearing in front 
of the integrals ( pT| ) disappears. We obtain an equation similar to the third line of Eq.(16), 
but this time containing a sum over four eigenvalues: 



m=l fc=l— 4 



e 



Wba - nX V {^^b2 + mX) {wbA + mX) (cvbi + mX) {wbz + mX) 
g,(,_,„^^,^^,^^^ (-p/A^)"e-"^'-° m^ + 2Ebm/X + Sb 



Wfea/A + n J^J-^ (wfei/A + m){wb2/\ + m){uJb3/>^ + m){wbA/X + m) 



TT {-xi + m){-X2 + m) 

(w6i/A + m){wb2/X + m){ujbz/X + m){wbA/\ + m) 

In the previous expressions we have used the following symbols: Eb = H^^bb- The factors 
Sb are respectively si^a = C3(/c2 — kf)/X'^ and S2,4 = C'i(/cf — fe^/A^. xi^2 are the roots of 



11 



the second order polynomial in m appearing in the third line; the roots xi^2 depend on b 
but not on a. 

To get an explicit expression for the diagonal elements of U, we take into account the 
definitions of xi, X2 and the quantities Wbb] in this way we obtain 



Ijin) ^ ^E,(r-ro) i-p/>^^Te (-X2)(n) ^^^^ 

(^^f.a/A)(n)(m/3/A)(„)(wfe^/A)(„) ■ 

The symbols a, (5,^ represent here the three non-repeated indices not equal to b. 
For the non-diagonal elements, we use the fact that for any b (for any row of the 
matrix) the following relation is true: 



XlX2 Wbb 

Hence, we have: 

Wba/A + n (n - 1)! (a;6«/A)(„)(u;b/3/A)(„)(a;b^/A)(„) ■ 

Once we have obtained Eqs.( ^j50| ), it is straightforward to perform the summation of 
the perturbative series in terms of generalized hypergeometric functions as was done for 
the one channel case. The evolution operator for the system given by Eq.(p^ is defined 
in terms of the auxiliary matrices Ured, Us as: 

U{r ^ oo,ro) = ex.p{HDr) Ured{ro), (52) 

with the matrix Ured given by : 

Ured{ro) = exp{-HDro) Usiro). 

The expressions for the elements of the matrix Us in terms of hypergeometric functions 
are easily deduced by inspection; for any row "6" the diagonal and non-diagonal elements 
are: 

Us,bb = 2F3{-Xbl,-Xb2;UJba/>^,^^bi3/>^,^b-y/^;z), (53) 
Us,bj{j^b) = Wbj/Wbb[Us,bl-2F3{-Xbl,-Xb2;l+UJbj/X,u;bl3,^b'y]z)] , (54) 

where coij = {H£))ii — {H£))jj and a, (3 and 7 represent non-repeated indices as before. 
Moreover Xbi,Xb2 are the roots of the equation: 

x'' + 2^x + Sb = 0, (55) 
A 

where the Sb are the quantities defined previously. 
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The evolution operator for a finite interval is the product of two matrices Ured defined 
at the first and last points of the interval: 

U{r,ro) = Ured{r)-^Ured{ro). (56) 

(2) 

In the "physical" basis, where Hq ' is diagonal, the general solution with arbitrary 
initial conditions at r = tq is given by: 

\ yi{r) J '^"'^ \ yi{ro) J 

By using the expression for w, a general solution with regular behavior at origin is given 
by: 

yiir) = (U[u] - U[i2] + Ui2i] - U[22]) w^u]R; (58) 

where R is an arbitrary constant two dimensional vector. We have denoted with U^^j^ 
(1 ^ hj ^ 2) the 2x2 block element of the matrix U. 

The asymptotic behavior for yi{r) is obtained from Eg. ([5^): 

yi{r — > oo) = (exp(zgr)e + exp{—iqr)a) {iq)^^R; (59) 

with the matrices: 

e = - Us,[12]) ; cr = {U,pi] - U,^[22]) ■ (60) 

According to Eq.(^), the scattering matrix is given in terms of the matrices e, o" as: 

Si{q) = ae-i . (61) 

The amplitude matrix K is obtained by expanding the exponentials appearing in 
Eq.(^) in trigonometric functions. Hence Eq.(|59D is equivalent to the following: 

yi{r — > oo) = {cos{qr){e + a) + i sin(gr)(e — a)) {iq)~^R. (62) 

Then, according to the usual definition, 

Ki = {€-a){e + a)-^. (63) 

It is easy to check that the standard relation between the scattering matrix Sj and the 
amplitude matrix Kj holds: 

Si = {l + iKi){l-iKi)~^. (64) 

The expressions for S and K in any other basis are obtained by a linear transformation. 
In the original basis where Eq.(^8|) was written, the matrices S^K are given by 

the matrix u^^^ is the 4x4 version of n, in a 2 x 2 block form is given by: 

""-(o !)■ («'^) 
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4 Conclusions 

In this work we have solved exactly an s-wave radial Schrodinger equation which in addition 
to a multichannel square-well potential contains a potential of exponential profile in one 
of the channels. The solution has been obtained by summing the infinite perturbative 
scries for the evolution operator of the system. The scattering matrix has been computed 
explicitly for the special case where the square-well and the exponential potential extend 
to infinity. The scattering length for the one channel case has been obtained in an explicit 
analytical form. Due to the factorizing property of the evolution operator, the results are 
easily applicable to situations in which short and long range potentials coexist. 
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